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Abstract 



QJ , We analyse the direct detection of neutralino dark matter in the framework 



X 



of the Next-to-Minimal Supersymmetric Standard Model. After performing a 
detailed analysis of the parameter space, taking into account all the available 
constraints from LEPII, we compute the neutralino-nucleon cross section, and 
compare the results with the sensitivity of detectors. We find that sizable val- 
ues for the detection cross section, within the reach of dark matter detectors, 
are attainable in this framework. For example, neutralino-proton cross sections 
compatible with the sensitivity of present experiments can be obtained due to 
the exchange of very light Higgses with m^o < 70 GeV. Such Higgses have a 
significant singlet composition, thus escaping detection and being in agreement 
with accelerator data. The lightest neutralino in these cases exhibits a large 
singlino-Higgsino composition, and a mass in the range 50 < m^o < 100 GeV. 



PACS: 12.60.Jv, 95.35. +d 



1 Introduction 



One of the most important enigmas in physics is the problem of the dark matter in 
the Universe. Particle physics, and in particular extensions of the standard model 
(SM) offer candidates for dark matter. Among the most interesting ones are Weakly 
Interacting Massive Particles ( WIMPs) , since these can be left over from the Big Bang 
in sufficient number to account for a significant fraction of the observed matter density. 

Since 1987, impressive experimental efforts have been carried out for the direct 
detection of WIMPs through elastic scattering with nuclei in a detector [1]. In fact, 
one of the experiments, the DAMA collaboration, reported data favouring the existence 
of a WIMP signal [2]. Once uncertainties on the halo model are taken into account 
[2, 3], this signal is compatible with WIMP masses smaller than 500-900 GeV and with 
WIMP-nucleon cross sections in the range u ^ 10~^ — 6 x 10~^ pb. 

However, this result has not been confirmed by the other collaborations. In partic- 
ular, CDMS Soudan [4], EDELWEISS [5] and ZEPLIN I [6] have excluded important 
regions of the DAMA parameter space^. In the light of these experimental results 
more than 20 experiments are running or in preparation around the world. For ex- 
ample, this is the case of GEDEON [8], which will be able to explore positively a 
WIMP-nucleon cross section cr > 3 x 10^^ pb. CDMS Soudan will be able to test in 
the future cr > 2 x 10~^ pb, and the very sensitive detector GENIUS [9], will be able 
to test a WIMP-nucleon cross section a ^ 10~^ pb. In fact, already planned detectors 
working with 1 tonne of Ge/Xe are expected to reach cross sections as low as 10~^° pb 
[10]. 

Given this situation, and assuming that the dark matter is a WIMP, it is necessary 
to analyse the theoretical predictions for the WIMP-nucleon cross section. Obviously, 
the answer to this question depends on the particular WIMP considered. The leading 
candidate in this class of particles is the lightest ncutralino, Xi^ which appears in 
supersymmetric (SUSY) extensions of the SM [11]. The cross section for the elastic 
scattering of Xi on nucleons has been examined exhaustively in the context of the 
Minimal Supersymmetric Standard Model (MSSM) [1] . In particular, there are regions 
of the parameter space of the MSSM where the neutrahno-nucleon cross section is 
compatible with the sensitivity of present (and future) dark matter detectors. 

However, it is well known that the MSSM faces a naturalness problem - the so- 
-'^For attempts to show that DAMA and these experiments might not be in conflict, see Ref. [7]. 
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called /i problem [12] - arising from the presence of a mass term for the Higgs fields 
in the superpotential, fiHiH2. The only natural values for the fi parameter are either 
zero or the Planck scale. The first is experimentally excluded since it leads to an 
unacceptable axion once the electroweak (EW) symmetry is broken, while the latter 
reintroduces the hierarchy problem. There exist explanations for an 0{Mw) value for 
the 11 term, although all in extended frameworks [12, 13]. 

The Next-to-Minimal Supersymmetric Standard Model (NMSSM) [14-25] provides 
an elegant solution to the fi problem of the MSSM via the introduction of a singlet 
superfield S. In the simplest form of the superpotential, which is scale invariant and 
contains the SH1H2 coupling, an effective /i term is generated when the scalar compo- 
nent of S acquires a vacuum expectation value (VEV) of order the SUSY breaking scale. 
This effective coupling is naturally of order the EW scale if the SUSY breaking scale is 
not too large compared with Mw- In fact, the NMSSM is the simplest supersymmetric 
extension of the standard model in which the EW scale exclusively originates from the 
SUSY breaking scale. Another appealing feature of the NMSSM is related to the "httle 
fine tuning problem" of the MSSM, or equivalently, the non-observation of a neutral 
CP-even Higgs boson at LEP II. As shown in [16], in the context of the NMSSM the 
latter problem becomes less severe. Although the symmetries of the NMSSM may give 
rise to the possibility of a cosmological domain wall problem [17], this can be avoided 
by the introduction of suitable non-normalisable operators [18] that do not generate 
dangerously large singlet tadpole diagrams [26]. These additional operators can be 
chosen small enough as not to alter the low energy phenomenology. 

In addition to the MSSM fields, the NMSSM contains an extra CP-even and CP-odd 
neutral Higgs bosons, as well as one additional neutralino. These new fields mix with 
the corresponding MSSM ones, giving rise to a richer and more complex phenomenology 
[15, 19-22]. A very fight neutralino may be present [19]. The upper bound on the mass 
of the lightest Higgs state is larger than in the MSSM [20]. Moreover, a very fight 
Higgs boson is not experimentally excluded [21,22]. All these properties may modify 
the results concerning the neutralino-nucleon cross section with respect to those of the 
MSSM. 

In fact, in comparison with the MSSM, there are only a few works in the fiterature 
studying the direct detection of the fightest neutrafino in the NMSSM [23,24], as 
well as its relic density [25] . Thus, given the recent experimental results concerning the 
detection of dark matter, and in view of the appealing theoretical and phenomenological 
properties of the NMSSM, it is important to carry out an up-to-date analysis of the 
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neutralino-nucleon cross section in this framework. Moreover, the recently pubhshed 
FORTRAN code NMHDECAY [21] allows a precise calculation of the particle spectrum 
in the NMSSM, as well as a complete check of all the available experimental constraints 
from LEP, thus enabling a thorough study of a in the allowed parameter space of the 
NMSSM. 

The outline of the paper is as follows: In Section 2 we introduce the model, dis- 
cussing in particular its Higgs potential, Higgs and neutralino mass matrices, and the 
parameter space. In Section 3 we examine the relevant effective Lagrangian describing 
the elastic xj'ii^cleon scattering and its associated cross section. Section 4 is devoted 
to the presentation of the results for the x'J'^ucleon cross section in the NMSSM, tak- 
ing into account the relevant constraints on the parameter space from accelerator data. 
Our conclusions are given in Section 5. 

2 Overview of the NMSSM 

In this Section, we review some important features of the NMSSM. In particular, we 
discuss the Higgs and neutralino sectors of the model, presenting the tree-level mass 
matrices and mixings which are relevant for our analysis. We also discuss the theoretical 
and experimental constraints, and how these are reflected in the parameter space. 

2.1 Higgs scalar potential 

In addition to the MSSM Yukawa couphngs for quarks and leptons, the NMSSM super- 
potential contains two additional terms involving the Higgs doublet superflelds, and 
H2, and the new superfield S, a singlet under the SM gauge group SU{?>)c x SU{2)l x 
U{1)y, 

W = eij (1; Hi Q'u + Yd Hi d + Y^Hi U e) - e^jX S H{Hi + \i^S\ (2.1) 

where we take Hi = {H^,H^), ifj = (i7^,if2), i,j are SU(2) indices, and ei2 = 1. 
In this model, the usual MSSM bilinear fi term is absent from the superpotential, and 
only dimensionless trilinear couplings are present in W . However, when the scalar 
component of S acquires a VEV, an effective interaction fj,HiH2 is generated, with 
fi = X{S). 

As mentioned in the Introduction, the superpotential in Eq. (2.1) is scale invariant. 
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and the EW scale will only appear through the soft SUSY breaking terms in >Csoft, 
which in our conventions is given by 

—J^soh — ''TT'Q Q* Q + iri^ u* u + m|, d* d + rn-^ L* L + e* e 
+ ml^ H* Hi + ml^ H;H2 + m| S*S 

+ eij Yu Hi Q'u + Ad Hi d + A, Y, H\ U e + H.c.) 
+ (^-e,,\A^SHiHi + ^KA,S^ + }l.c^ 

- ]- {Ms A3 A3 + M2 A2 A2 + Ml Ai Ai + H.c.) . (2.2) 

In our subsequent analysis we assume that the soft breaking parameters are free at the 
EW scale. In addition to terms from £soft, the tree-level scalar Higgs potential receives 
the usual D and F term contributions: 

Vf - |Ap (li/iH^p + \H2\^\S\^ + \eijHlHi\^) + \f,\^\S\^ 

- {eijXK*HiHiS*^ + H.c.) . (2.3) 

2.2 Minimization of the tree level scalar potential 

Once the EW symmetry is spontaneously broken, the neutral Higgs scalars develop the 
following VEV's: 

{H',)^v,, {Hl)^V2, {S)^s. (2.4) 

One has to ensure the absence of non-vanishing VEV's for the charged Higgs fields, 
which would induce the appearance of charge breaking minima. By means of an 
SU{2)l X U{1)y transformation, one can take, without loss of generality, {H2) — 
while {H2) — V2 E M+. The condition to have Vi — (H^) = as a global minimum is 
quite involved; still, imposing that f ]~ = is a local minimum is equivalent to requiring 
that the physical charged Higgses have positive mass squared. 

Bringing together all the terms in Eqs. (2.2) and (2.3), we then obtain for the 
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tree-level neutral Higgs potential: 

I I \ |2 /l |2| |2 , I |2| |2 , I |2| |2\ , I |2| |4 

+ |A| [\s\ \vi\ + \s\ \v2\ +\vi\ |f2| j + \s\ 
+ "^i^il^iT + "^kl^2|^ + ml\s\'^ 

+ (^-Xk*viV2S*'^ - XAxsviV2 + ^kA^s^ + H.c.^ . (2.5) 

In the following, we assume that A, k, as well as the soft SUSY breaking terms are real. 
This implies the absence of explicit CP violation in the scalar sector. Although Vi and 
s can be complex parameters, the global Z3 symmetry exhibited by the supcrpotential 
implies that CP- violating extrema of V^^utfai maxima rather than minima [27]. In 
principle. A, k, and the trilinear soft-breaking terms. Ax and A,^, in Eq. (2.5) can have 
both signs. 

Ensuring that the tree-level potential has a minimum with respect to the phases 
of the VEV's directly excludes some combinations of signs for the parameters. After 
conducting this analysis, and given that the potential is invariant under the symmetries 
A, K, s — —A, — K, —s and A, vi —A, —Vi, we adopt, without loss of generality, the 
sign convention where both A and Vi are positive. We then have only positive values 
for A and tan/3, while k and jj, (= As), as well as Ax and A,^, can have both signs. 

In what follows, we summarise the conditions for k, A\. A^ and /i (= As) obtained 
from the minimization of the potential with respect to the phases of the VEV's. In 
particular, for k > 0, one can analytically show that minima of V^^utfai t>e ob- 
tained for the following three combinations of signs, provided that in each case the 
corresponding conditions are fulfilled, 

(i) sign(s) = sign(AA) = -sig-ii{A^), 

which always leads to a minimum with respect to the phases. 

(ii) sign(s) = -sign(^;^) = -sign(A«), 

with I^kI > 3AviV2|A;y|/(— |syl;y| + fi;|s^|), where the denominator has to be posi- 
tive. 

(iii) sign(s) = sign(^A) = sign(AK), 

with \A^\ < Z\viV2\Ax\/{\sAx\ + k\s^\). 

Similarly, for k < 0, minima can only be obtained for the combination 
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(iv) sign(s) = sign(AA) = sign{A^), 

with \ Ak\ > 3XviV2\Ax\/{\sAx\ — ac|s^|), where the denominator has to be positive. 

Numerically, one finds that these tree-level conditions generally hold even after the 
inclusion of higher order corrections. 

One can derive three minimization conditions for the Higgs VEV's and use them 
to re-express the soft breaking Higgs masses in terms of A, k, Ax, A^, vi, V2 and s: 

m]j^ = - A^ [s^ + v'^ sin^ /3) - -M| cos 2/3 + As tan /3 {ks + Ax) , 

= - A^ (s^ + cos^ 0) cos 2/3 As cot /3 (ks ^a) , 

m| = - )kv^ - 2k^s^ + Xkv"^ sin 2/3 + ^^^^ sin 2/3 - /«A«s , (2.6) 

where v'^ = vl-\-vl = 2M'^/gl and tan/3 = ^2/^1. 



2.3 Higgs boson mass matrices 

Subsequent to EW symmetry breaking, and after rotating away the CP-odd would-be 
Goldstone boson, we are left with five neutral Higgs states and two charged Higgs 
states. Assuming 

rrO _ , + iHii _ H2R + iH2i Sr + iSi 

^ " ^ 7! ' ^2 = ^2 + ^ , S = s + — -^=— , (2.7) 

among the neutral Higgses we find three CP-even states - Hir, H2R, Sr and two CP-odd 
components, 5"/, with A^ related to the original fields as i?i(2)/ = sin /3(cos /3) Us- 
ing the minimization conditions above, the tree-level mass matrix for the neutral Higgs 
bosons can be easily obtained. Since we have made the assumption that there is no 
CP-violation on the Higgs sector, CP-cvcn and CP-odd states do not mix, and the cor- 
responding mass matrices can be written in the respective basis, — {Hir, H2R, Sr) 
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and P*^ = {A^, Sj). For the CP-even states, we have 

M%ii = M| cos^ (3 + Xs tan (3(Ax + ks) 
Al| 22 = sin^ /3 + As cot ^(Ax + ks) 

s 

( \ 
M\^2 = ( A^^^ - ^ J sin 2/3 - As {Ax + fi;s) 

A^l 13 = 2A^vis - Xv2 {Ax + 2«;s) 

Ml^22. = 2A^V2S - A^;i {Ax + 2ks) . (2.8) 

The CP-even Higgs interaction and physical eigenstates are related by the transforma- 
tion 

K = Sa,H',, (2.9) 

where 5* is the unitary matrix that diagonalises the above symmetric mass matrix, 
a,b — 1,2,3, and the physical eigenstates are ordered as^ m^o < rUf^o < mf^o. In 
the pseudoscalar sector, after rewriting the CP-odd mass terms in the P° basis, the 
corresponding (symmetric) mass matrix reads 

M%22 ^x(^K + sin 2/3 - 3kA^s 

M%i2 = Xv {Ax - 2ks) , (2.10) 
and the relation between physical and interaction eigenstates is given by 

<-PijPl (2.11) 
Regarding the charged Higgs mass, at the tree level it is given by 

m^± =M^-AV + A(AA + fi:s)^^. (2.12) 

2.4 Neutralino mass matrix 

When compared to the MSSM case, the structure of chargino and squark mass terms 
is essentially unaffected, provided that one uses fi = Xs. However, in the neutralino 
^Throughout the paper we always adopt the convention rrij < m,- for i < j. 
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sector, the situation is more involved, since the fermionic component of S mixes with the 
neutral Higgsinos, giving rise to a fifth neutralino state. In the weak interaction basis 
defined by \E'°^ = (^B^ = —iX', = —iXs, H^, H2, , the neutralino mass terms in 
the Lagrangian are 



(2.13) 



with M.ySi a 5 X 5 matrix, 



Ml 






M2 



— M^ sin sin /? \ 
Mz cos dw cos /? — Mz cos Qw sin /3 



V 



— M2 sin^T^ cos/3 cos ^^vF cos /3 

M^ sin sin /? — M^cos^vK sin/3 





-As 



-As 


— Awi 



The above matrix can be diagonalised by means of a unitary matrix N ^ 

N* M.^oN~^ = diag(m^o, m^o, m^o, m^o, m^o) , 



— AW2 

— Aiii 
2ks I 
(2.14) 



(2.15) 



where m^o is the lightest neutralino mass. Under the above assumptions, the lightest 
neutralino can be expressed as the combination 



(2.16) 



In the following, neutralinos with N^^^ + Nl/^ > 0.9, or Nf^ > 0.9, will be referred to as 
Higgsino- or singlino-like, respectively. 



2.5 NMSSM parameter space 

At the weak scale, the free parameters in the Higgs sector are (at tree level): A, k, 
m^fj^, rrijj^, m|. Ax and A,^. Using the three minimization conditions of the Higgs 
potential (including the dominant one- and two-loop corrections), one can eliminate 
the soft Higgs masses in favour of M^,tan/5 and 11. We thus consider as independent 
parameters the following set of variables 

A, tan/3, ^A, (2.17) 

In our study, the soft scalar masses as well as the soft gaugino soft masses Mj are 
free parameters at the EW scale. We scanned over the parameter space using the 
program NMHDECAY [21] and in what follows we overview the most relevant aspects of 
the analysis. 
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For each point in the parameter space, one requires the absence of Landau singular- 
ities for A, K, Yf and Yf, below the GUT scale. For m^"'*^ = 175 GeV, this translates into 
A < 0.75, \k,\ < 0.65, and 1.7 < tan/? < 54. In addition one verifies that the physical 
minimum is a true one, in other words, that it is deeper than the local unphysical 
minima with {H'(^ = and/or (S) = 0. 

One then computes the scalar, pseudo-scalar and charged Higgs masses and mix- 
ings, taking into account 1- and 2-loop radiative corrections. The dominant 1-loop 
corrections to the Higgs masses originate from top, stop, bottom and sbottom loops, 
and the corresponding corrections to m^,, are of 0{Y^\). Pure electroweak contribu- 
tions of 0{g^) are also taken into account. Regarding 2-loop corrections to the effective 
potential, the dominant ones are associated with top-stop loops, and only the leading 
(double) logarithms are included. The chargino and neutralino masses and mixings 
are computed and the couplings of the scalar and pseudoscalar Higgs to charginos and 
neutralinos are calculated. 

Finally, all available experimental constraints from LEP are checked: 

1) In the neutralino sector, we check that the lightest neutralino does not contribute 
excessively to the invisible width of the Z boson {T(Z — >• XiXi) < 1-76 MeV [28]) if 
m^o < Mz/2, and that a{e~^e~ XiXi) < 10"^ pb if m^o + m^o < 209 GeV (i > 1) 
and a{e+e- x°X°) < 10"^ pb if m^o + m^o < 209 GeV > 1) [29]. 

2) In the chargino sector, we verify that the lightest chargino is not too fight {m^+ > 
103.5 GeV [30]). This leads to a lower bound on |//| > 100 GeV. 

3) In the charged Higgs sector, we impose mH+ > 78.6 GeV [31]. 

4) In the neutral Higgs sector, we check the constraints on the production rates (reduced 
couplings) X brancfiing ratios versus the masses, for all the CP-even states /i'^ and CP- 
odd states a°, in all the channels studied at LEP [32]: e+e" h?Z, independent of the 
hP decay mode (IHDM); e+e" — > hPZ, dependent on the hP decay mode (DHDM), with 
the Higgs decaying via hP hh, hP hP 2 jets hP 77 and hP invisible; 
associated production modes (APM), e+e~ — > hPa^, with hPaP — > 4fe's, hPa^ 4t's 
and hPaP — > a°a°a° — > 66's (see [21] for a detailed discussion). 

It is worth noticing here that other available experimental bounds such as e.g. the 
6 — > 57 branching ratio, and the current upper limit on tfie decay Bg — > , migfit 
also produce relevant constraints on tfie parameter space. Tfie presence of ligfit Higgs 
states, in particular tfiat of very ligfit pseudoscalars, migfit translate into potentially 
large contributions to tfiese processes. Tfiis issue fias begun to be addressed in [33] 
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for the large tan (3 regime. For instance, the branching ratio Bs ji'^jj,' scales as 
tan^/3/m^o, thus producing important constraints [34] in non-universal MSSM sce- 
narios associated with light pseudoscalar Higgs, which favour large neutralino-nucleon 
cross sections. These constraints become specially relevant in the large tan/3 regime 
(tan/? > 35). However, notice, that although in the NMSSM scalar and pseudoscalar 
Higgs can indeed be very light, large neutralino-nucleon cross sections are in general 
obtained for the low tan /5 regime, as we will discuss in the next sections. It is impor- 
tant then to carry out a detailed analysis of these processes. In particular, 6, and 
B decays, together with the muon anomalous magnetic moment, may play an impor- 
tant role in further constraining the NMSSM parameter space. Such a study will be 
considered in a forthcoming publication [35] . 



3 Neutralino-nucleon cross section 

The most general supersymmetric low-energy effective four-fermion Lagrangian that 
describes the elastic scattering of the lightest neutralino with the nucleon is given 
by [36,37] 

^eff =x\ Yl5 Xl Qi 7m ("li + "2i75) Qi + "Si Xl Xl Qi Qi 

+ «4i Xl 75 Xl Qi 75 Qi + «5i Xl Xl Qi 75 Qi + "ei X? 75 Xi Qi Qi: (3-18) 

where i — 1,2 denotes up- and down-type quarks, and the Lagrangian is summed over 
the three quark generations. In the absence of CP-violating phases, the terms propor- 
tional to as and vanish. Moreover, those associated with cti and (as well as 
and cte, should these be present) are velocity-dependent, and can be safely neglected for 
our present purposes. The cross section associated with the spin-dependent coefficient 
(q;2) is only non-zero if the target nucleus has a non-vanishing spin and, contrary to 
case of the scalar (spin-independent) term, adds incoherently. For the case of heavy 
targets, as those used in the experiments mentioned in the Introduction, the scalar 
cross section associated with is in general substantially larger, and henceforth we 
shall focus our discussion on the latter. 

We begin by decomposing a^i into two terms, one arising from squark s-channel ex- 
change and the other from the t-channel, neutral Higgs mediated interaction (Fig. 1 
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(a) (b) 

Figure 1: Feynman diagrams contributing to the neutralino-nucleon scalar cross section: 
(a) squark exchange and (b) scalar Higgs exchange. 

(a) and (b), respectively). One obtains^ 

4 = -i: 4(„^.L.„) R-[(g.")(Cf)1. (3.19) 

X=l ^ XV 

c^^ = J2^^y^''^^HL] , (3.20) 

III' L n 
a=l hO 

where X = 1,2 denotes the squark mass eigenstates and a = 1,2,3 refers to the 
scalar Higgs mass eigenstates. The relevant NMSSM couplings for the neutralino- 
squark-quark (C^}j), neutralino-neutralino-Higgs (C^^_^) and Higgs-quark-quark (Cy) 
interactions are given in Appendix A. 

The term is formally identical to the MSSM case, differing only in the new 
neutralino mixings stemming from the presence of a fifth component. In particular, 
in regions of the NMSSM parameter space where the singlino component dominates 
the lightest neutralino state, there will be a significant reduction in the Bino- and 
Wino-squark-quark couplings, and hence in al^. 

Regarding the Higgs mediated interaction term (a^i), the situation is slightly more 

involved since both vertices and the exchanged Higgs scalar significantly reflect the 

new features of the NMSSM. First, let us recall that in regions of the parameter space 

where the lightest Higgs boson has a sizable singlet component , the Higgs-quark-quark 

coupling might be substantially reduced. Regarding the Higgs-Xi-Xi interaction, in 

■^When compared to the analogous expression of Rcf. [24], we find some discrepancies in a^^, namely 
a missing singlet- Higgsino-Higgsino term (proportional to A), and the presence of an additional overall 
weak coupling constant, g. 
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addition to a new component in the lightest state, the most important alteration 
emerges from the presence of new terms, proportional to A and n (cf. Appendix A). 
Nevertheless, and as already mentioned, light Higgs bosons can be experimentally 
allowed in the context of the NMSSM. Should this occur, and if these states are not 
pure singlets (thus displaying a non-vanishing coupling to matter) the exchange of 
light Higgs scalars in the t-channel might provide a considerable enhancement to the 
neutralino-nucleon cross section. 

It is worth mentioning that an enhancement of a^^ with respect to a'^- is not an 
effect unique to the NMSSM. In fact, it has been already noticed that in the MSSM, 
and once the mSUGRA inspired universality for the soft scalar and gaugino masses 
is abandoned, the cross section associated with the channels involving scalar Higgs 
exchange can be substantially enhanced. Similar to what will occur in the present 
model, the MSSM t-channel contributions become larger once the Higgsino components 
of x5 are augmented and/or the Higgs masses are reduced [39] (e.g. via non-universal 
soft masses at the GUT scale [40]). 

The scalar interaction term contributes to the x5"Nucleon cross section as 

asN ^^f^, (3.21) 

TT 

where is the Nucleon-x5 reduced mass, = mNm^o/{mN + ?7t.^o), and 

q=u,d,s ^ q=c,b,t ^ 

In the above, rriq is the quark mass, and the parameters f!j^^ are defined as {N\mqqq\N) 
= m^f?^ /(^^ can be derived from ff^^ as /^^ = 1 - E,=m,sC- Fallowing [38], 
we take the following values for the hadronic matrix elements: 

/jfj = 0.020 ± 0.004 , = 0.026 ± 0.005 , f^^] = 0.118 ± 0.062 , 

/j,^) = 0.014 ± 0.003 , f^J = 0.036 ± 0.008 , f^J = 0.118 ± 0.062 . (3.23) 

In the numerical analysis of the next section we will use the central values of the above 
matrix elements. Notice that f^J = f!j^'^ and both are much larger than J'tq for u 
and d quarks, and therefore fp and /„ are basically equal. Thus we will focus on the 
neutralino-proton cross section, 

CTsp = cr^o-p = fp , (3.24) 



with TJir = mpTJi^o I (mp + m^o) ~ rrip. 
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4 Results and discussion 



In this Section the viabihty of the detection of the hghtest NMSSM neutrahno as a 
dark matter candidate, will be studied. In particular, we will compute the theoretical 
predictions for the direct detection of neutrahnos through their elastic scattering with 
nucleons inside a detector. In our computation we will take into account relevant 
constraints on the parameter space from accelerator data. On the other hand, given the 
complexity of the computation of the relic neutralino density, we prefer to consider in a 
forthcoming publication [41] the constraints arising from reproducing the WMAP data 
[42], 0.094 < floyih^ ^ 0.129, on our relevant parameter space for the cross section. In 
particular, we will see in this section that light pseudoscalars are present in interesting 
regions of the parameter space, and this might translate into large contributions to 
the annihilation cross section, implying a reduction in the associated relic density. 
Therefore we may expect the WMAP lower bound to play an important role in those 
points with a very large neutralino-nucleon cross section. For this study we will use 
the general analysis of the relic neutralino density, including coannihilations, that has 
been carried out for the NMSSM in [43]. 

As discussed in the Introduction, many experiments for the direct detection of dark 
matter are running or in preparation. Thus, in our analysis, we will be particularly 
interested in the various NMSSM scenarios which might potentially lead to values of 
cr^o_p in the sensitivity range of those detectors. 

Although the free parameters in our model have already been presented in Sec- 
tion 2.5, it is worth recalling that the Higgs and neutralino sectors of the theory are 
specified by 

A, K, //(=As), tan/3, A^, A^, Mi, M2 . (4.25) 

As aforementioned, we take these parameters to be free at the EW scale. Based on an 
argument of simplicity"', the low-energy squark masses and trilinear couplings, which 
appear in the computation of the neutralino-nucleon cross section, are taken to be 
degenerate^. Unless otherwise stated, the common SUSY scale will be Msusy = 1 TeV. 

^Since in our analysis of the neutralino-nucleon cross section the detection channels mediated by 
Higgs scalars will be enhanced with respect to those mediated by squarks, the sensitivity of the results 
to variations of the squark parameters will be very small. 

■''Regarding the stop mass matrix we will work in the maximal-mixing regime, where the off-diagonal 
term takes the form nitXt = mt a/S -Msusy- Departures from this case would not affect significantly 
the theoretical predictions for the neutralino-nucleon cross section. 
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Having free squark and slepton soft parameters at the EW scale allows us to ensure 
that in our analysis the lightest SUSY particle is indeed the Xv Also led by arguments 
of naturalness, we shall take a lower bound for A, \„im ~ /u/smax- Thus, taking the 
conservative range s ^ 10 TeV, this translates into Amin ~ /x(GeV) x 10"'^. 

We begin our analysis by taking values for the soft gaugino masses that mimic at low 
scale the results from a hypothetical unified value at the GUT scale. Consequently, we 
will choose M2 = 1 TeV and Mi = 500 GeV. For the gluino mass, the value M3 = 3 TeV 
will be taken. Later on we will address variations of these values. In the following we 
take \n\ > 110 GeV, since in most cases this allows to safely avoid the LEP bound on 
the lightest chargino mass. Throughout this Section, we shall consider several choices 
for the values of A\, A^;, /i and tan/5, and for each case, we will study the associated 
phenomenology. 

In order to simultaneously analyse the dark matter predictions and understand the 
effect of the experimental constraints on the parameter space of the NMSSM, it is very 
illustrative to begin our study in the plane generated by the Higgs couplings in the 
supcrpotential, A and k,. In Section 2.2 we commented on the conditions to be applied 
to each of the sign combinations of the parameters, which arise from ensuring that the 
tree-level potential has a minimum with respect to the phases of the VEV's. In the 
following it will be clarifying to discuss each case separately. Let us first consider the 
cases associated with positive values of k. 

4.1 iiA,^ < and /lAx > > 0) 

As a first choice, we will consider the two cases where < and iJ,Ax > 0, namely 
those with /i, Ax, — > and /x, ^4;^, —A^ < 0. 

In both cases, part of the parameter space can be excluded due to the occurrence 
of tachyons in the CP-even Higgs sector. Namely, it is easy to see from the expression 
of the CP-even Higgs matrix (2.8) that the off-diagonal terms |A^|i3| or |A^|23l '^^^ 
become significantly bigger than TWigg, thus leading to the appearance of a negative 
eigenvalue. This will typically happen for moderate to large values of A and small k,, for 
which m^o is small. Large values of \Af^\ and tan (3 lead to an increase of the tachyonic 
region, as we will later see. On the other hand, the eigenvalues of the CP-odd Higgs 
mass matrix are never negative. The CP-odd Higgs masses also decrease for large A 
and small k, but their minimum value is bounded by the appearance of tachyons in 
the CP-even sector. 
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Figure 2: {X,k) parameter space for tan/3 — 3, A\ — 200 GeV, — —50 GeV and 
fjL — 110 GeV. In both cases, the ruled area represents points which are excluded due to 
the occurrence of a Landau pole. The grided area is excluded because of the appearance of 
tachyons. The grey area is associated to those points that do not satisfy the LEP constraints 
or where (at least) the LEP bound on direct neutralino production is violated. Dotted lines 
in the experimentally accepted region represent contours of scalar neutralino-proton cross 
section cr^o_p. In (a), from top to bottom, solid lines indicate different values of lightest 
Higgs scalar mass, m^jO = 114, 75, 25 GeV, and dashed lines separate the regions where 
the lightest scalar Higgs has a singlet composition given by 5*^1 = 0.1, 0.9. In (b), from 
top to bottom, solid lines are associated with different values of the lightest neutralino 
mass, m^o = 100, 75, 50 GeV, while dot-dashed lines reflect the singlino composition of 
the lightest neutralino, iVfg = 0.1, 0.9. 

The (A, k) parameter space is shown in Fig. 2 for an example with tan/9 = 3, 
Ax = 200 GeV, = —50 GeV and ^ = 110 GeV. The points which are excluded 
due to the occurrence of a Landau pole are indicated, as well as those not fulfilling the 
experimental constraints. According to the discussion above, the tachyons appearing 
in the lower right corner are due to the CP-even Higgs sector. It is worth remarking 
that in these cases, due to the smallness of the lightest CP-even Higgs mass, the 
experimental constraints (see Section 2.5) from e+e~ h'^Z, both IHDM and DHDM 
(/i° — > bb, hP — > T+T~, and hP — > 2 jets), become very important and typically exclude 
the regions in the vicinity of those excluded by tachyons. 
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Dashed lines in Fig. 2a indicate the singlet composition of the lightest scalar Higgs. 
Singlet-like Higgses can be found for small values of /t, whereas doublet-like Higgses 
appear for large k. This can be qualitatively understood from the expression of the 
corresponding mass matrix (2.14). In particular, the diagonal term ^^533 becomes 
very small when n decreases. Interestingly, when the singlet composition is significant, 
the reduced coupling can be smaller and thus Higgses with m/jO < 114 GeV can escape 
detection and be in agreement with experimental data. This opens a new window in 
the allowed parameter space, characteristic of the NMSSM, which can have relevant 
consequences for dark matter detection as we will discuss below. 

In Fig. 2b the same case is represented, but emphasizing the information on the 
neutralino properties. The singlino composition of the lightest neutralino is shown 
with dot-dashed lines, while solid lines correspond to different values of its mass. As 
one would expect from the structure of the neutralino mass matrix (2.14), for small 
K, the lightest neutralino is essentially a singlino, with a small mass which can be 
approximated as m^o ~ 2^k/\. In the present case, singlino-like neutralinos appear 
for K < 0.04 and A < 0.2, whereas heavier, Higgsino-like, neutralinos (due to our choice 
of input values with /i < Mi) populate the rest of the parameter space. Regions with 
small masses of the neutralino may be excluded due to the bound on direct neutralino 
production, which becomes quite severe for light Higgsino-like neutralinos. 

In both figures, the different values of the neutralino-nucleon cross section are repre- 
sented with dotted lines. As already commented in Section 3, the cross section increases 
in those regions with a light CP-even Higgs, as long as it is not a pure singlet. This 
behaviour is clearly illustrated in these figures, which feature very large values for cr^o_p 
in the vicinity of the areas where the hghtest Higgs becomes tachyonic. On the other 
hand far from these regions the cross section stabilizes at 10~^ pb > cr^o_p > 10~^ pb. 

In order to illustrate this point in more detail, we have represented in Fig. 3 the 
resulting o"^o„p versus the lightest Higgs mass and the neutralino mass. Black dots fulfil 
all the experimental constraints, whereas grey dots arc those experimentally excluded 
(we do not plot those regions ruled out due to theoretical arguments, such as the 
occurrence of a Landau pole). The sensitivities of present and projected dark matter 
experiments are also depicted as a function of m^o for comparison. The small area 
bounded by dotted fines is allowed by the DAMA experiment in the simple case of 
an isothermal spherical halo model. The larger area also bounded by dotted fines 
represents the DAMA region when uncertainties to this simple model are taken into 
account. For the other experiments in the figure only the spherical halo has been 
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Figure 3: Scatter plot of the scalar neutralino-nucleon cross section, (T^o_p, as a function 
of (a) the neutralino mass, m^o, and (b) the lightest scalar Higgs mass, m/jO, for ^4;^ = 200 
GeV, ijl — 110 GeV, — —50 GeV, and tan/3 = 3. Black dots correspond to points 
fulfilling all the experimental constraints, whereas grey dots represent those excluded. In 
(a) the sensitivities of present and projected experiments are also depicted with solid and 
dashed lines, respectively. The large (small) area bounded by dotted lines is allowed by the 
DAMA experiment when astrophysical uncertainties are (are not) taken into account. 

considered in their analyses. In particular, the (upper) areas bounded by solid lines 
are excluded by EDELWEISS^ and CDMS Soudan. Finally, the dashed lines represent 
the sensitivities of the projected GEDEON, CDMS Soudan, and GENIUS experiments. 

Very large values for the cross section could in principle be obtained. However, 
these are associated to very light Higgses and are therefore subject to the strong con- 
straints on e+e" —>■ hPZ discussed above. Once every constraint is taken into account, 
points with cr^o_p < 10~^ pb appear, which correspond to light scalar Higgses with 
m^o > 75 GeV, surviving the experimental constraints due to their important singlet 
character, S^^ > 0.85. This is a clear consequence of the NMSSM that we will exploit in 
subsequent examples, since it allows for a significant increase in the cross section. The 
neutralinos in these regions have m^o > 70 GeV and have a mixed singlino-Higgsino 

^ Since the exclusion area due to ZEPLIN I is similar to EDELWEISS we have not depicted it here, 
nor in any subsequent Figure. 



18 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 50 100 150 200 250 

A m, (GeV) 



Figure 4: (A, k) parameter space with the corresponding constraints and neutralino-nucleon 
cross section as a function of the lightest neutralino mass for the case Ax — 200 GeV, 
— —200 GeV, fx — 110 GeV and tan/3 = 3. In the (A, k) plane the mass and 
composition of the lightest scalar Higgs, the composition of the lightest neutralino (only 
the line with N^^ = 0.1), and the predictions for cr^o_p are represented with the same 
line conventions as in Fig. 2, and the new ruled area (vertical lines) is excluded due to the 
occurrence of unphysical minima. The colour convention for the plot cr^o_p versus m^o is 
as in Fig. 3. 

composition (A^fg < 0.3 and N^^ + -^14 ^ 0.7 in the region with larger cross section). 

Notice that cr^o_p displays an important suppression around m^o 90 GeV. This is 
due to the cancellation of the contribution of the cross section coming from neutralino- 
neutralino-Higgs interaction due to the occurrence of terms with different signs. This 
type of accidental cancellations is analogous to those appearing in MSSM analyses for 
/X < [38]. 

These results are also sensitive to variations in the rest of the input parameters (A^, 
tan/3, /I, and Ax). For instance, increasing (i.e., making it more negative) leads 
to a further decrease in A4s,33 CP-even Higgs mass matrix, and therefore lighter 

Higgses can be obtained with a larger singlet composition. Although this translates into 
an enlargement of the regions where one has a tachyonic scalar Higgs, one may never- 
theless find a larger cr^o_p in the allowed areas. Choosing tan P — 3, Ax — 200 GeV and 
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/X = 110 GeV, but with = —200 GeV, one can obtain cr-^o^p > 10^"^ pb (points which 
in fact are already excluded by direct dark matter searches). The corresponding (A, k) 
parameter space, as well as o-^o^p versus the neutralino mass, are represented in Fig 4. 
Remarkably, very light Higgses are allowed in this case (m^jO > 20 GeV) due to their 
significant singlet character (0.9 < S^^ < 0.95). Once again, the hghtest neutralino 
exhibits a large singhno-Higgsino composition (iVfg < 0.3 and N'}^ + N^^ > 0.7). For 
these reasons, one hardly finds experimentally excluded regions: only narrow stripes, 
mostly due to direct production of x° and hP — > hh. Also, for small values of A and k, 
a very thin region excluded by the existence of false minima (see Section 2.5) appears. 
Conversely, decreasing \Ak\ helps reducing tachyonic regions. In the particular case 
where A,^ — 0, no tachyons emerge from the CP-even sector. The implications of this 
variation in the value of cr-fO_„ are minimal. 

Changing tan {3 has an important impact in the analysis, mainly due to the effect 
on the Higgs sector. The tachyonic regions become larger as tan (3 increases (extending 
towards higher values of A and k). As a consequence, the neutralino is never a pure 
singlino and its mass increases due to the larger mixing with Higgsinos. For this reason 
the exclusion due to direct neutralino production becomes larger. In the end, not only 
the allowed region is reduced, but also the predictions for (J^o_p are smaller. Also, for 
very small values of tan/3 very light Higgsino-like neutralinos can be found in large 
regions of the parameter space. The experimental constraints are, nevertheless, more 
important and only small areas survive. Fig. 5 illustrates these properties. 
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Figure 5: The same as in Fig. 4 but for the cases Ax — 200 GeV, A,^ — —50 GeV, // = 110 
GeV, and tan/3 = 2, 5, 10, from top to bottom. In the case with tan/3 = 2, only the lines 
with m/jO = 75, 25 GeV are represented, since m/jO < 110 GeV. 




Figure 6: The same as in Fig. 3a but for the cases A\ = 200 GeV, — —200 GeV, 
fj, — 110 GeV, and tan/? — 2, 4, 5, from left to right and top to bottom. 

Similar examples, but for = —200 GeV can be found in Fig. 6, where the predic- 
tions for o"^o_p are depicted as a function of the neutralino mass for tan/? = 2, 4, 5. As 
already mentioned, small values of tan (3 favour lighter neutralinos with larger detection 
cross section. 

Heavier neutralinos with a larger singlino composition can be obtained if the value 
of II is increased. For this reason, the regions where direct neutralino production is not 
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Figure 7: The same as in Fig. 3a but for the cases Ax — 200 GeV, tan j3 — 3, /j, — 200 GeV, 
and A« = -50, -200 GeV, from left to right. 

in agreement with experimental bounds become much narrower. The mass of scalar 
HiKKses also grows in this well as their doublet character. Constraints on 

the Higgs sector are still strong enough to forbid those points where the neutralino is 
mostly singlino, and in the end Xi preserves its mixed singlino-Higgsino character. In 
the remaining allowed area the predictions for o'^o_p can vary, being typically smaller 
than in cases with low /j. This is shown in Fig. 7 for two examples with fi = 200 GeV, 
Ax = 200 GeV, tan/? = 3, and = —50, —200 GeV. In particular, in the case with 
= —200 GeV, the detection cross section is much smaller than in the analogous 
example with n = 110 GeV presented in Fig. 4. 

Finally, variations in the value of Ax also influence the theoretical predictions on 
Cx°-p- There is a range of Ax for which the eigenvalues of the CP-even Higgs mass 
matrix are positive. However, for smaller or larger Ax, off-diagonal terms may become 
large enough to ease the appearance of tachyons in the large A regime. For instance, in 
Fig. 8 we have represented the (A, k) plane and the corresponding predictions for cr^^o-p 
in two cases with A^ = —50 GeV, fi = 110 GeV, tan/5 = 3, and Ax = 50, 450 GeV. 
We find that, in agreement with the discussion above, the tachyonic regions are larger 
than those for Ax = 200 GeV in both cases. Also, the areas excluded by experimental 
constraints associated to IHDM and DHDM are more extensive, and in the case of 
Ax = 450 GeV they forbid most of the parameter space. The neutralino is mostly 
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Figure 8: The same as in Fig. 4 but for the ca 
and Ax — 50,450 GeV, from top to bottom. 
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Higgsino in the remaining allowed points, with Nf^^ < 0.1(0.2) and A'^^3 + A'^^4 > 0.9(0.8) 
for Ax = 50(450) GeV, and there is a slight decrease in the predictions for cr^o_p. 

The range of values of Ax for which the allowed area is more extensive is very 
dependent on the rest of the inputs. In particular, since large tan/3 and \Ak\ increase 
the diagonal term, M.'j-,^^. in the CP-odd Higgs mass matrix, larger values of Ax can 
be taken before A1pi2 g^ts too big. For example, in the case with A^ — —200 GeV one 
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Figure 9: The same as in Fig. 3a but for the cases = —200 GeV, ji — 110 GeV, 
tan/3 = 4 and Ax = 300 GeV, on the left, and = -200 GeV, = 110 GeV, tan/3 = 5 
and = 450 GeV on the right. 

can still obtain large accepted regions for A\ = 300 — 450 GeV and tan /3 = 4 — 5, as 
evidenced in Fig. 9, where points entering the sensitivities of the present dark matter 
detectors are obtained with m^o < 75 GeV. 

To complete the analysis of the cases with /lAx > and < 0, we must address 
the possibility of having //, Ax, —A^ < 0. Note from (2.5) that the tree-level potential, 
K^utoi' therefore the Higgs mass matrices, are invariant under the exchange of 
the signs of fi, Ax and A^, provided that the signs of fxAx and /xA^ do not change. 
This implies that the above analysis regarding the Higgs sector is identical in this case. 
Differences arise, however, in the neutralino sector since the signs of Mi^2 were not al- 
tered. Therefore, the neutralino mass spectrum differs, as well as the lightest neutralino 
composition. Also the experimental constraints exhibit a slight variation. This case 
presents the same qualitative behaviour as the one formerly discussed in what the min- 
imization of the Higgs potential is concerned. Nevertheless, differences arise regarding 
the theoretical predictions for o'^Q_p due to the experimental constraints and the differ- 
ent position of the accidental suppressions in the Higgs-exchangc diagrams. These dif- 
ferences can be sizable for large tan/3. For instance, we have represented in Fig. 10 two 
examples with Ax = -200 GeV, = -110 GeV, tan/3 = 3, and A^ = 50, 200 GeV, 
where the suppression in cr^o_p is found to occur for m^o pa 110 GeV. 
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Figure 10: The same as in Fig. 3a but for the cases Ax — —200 GeV, /j, — —110 GeV, 
tan/3 = 3, and = 50, 200 GeV, from left to right. 

To sum up, we have found that large values of CTj^o-p, even within the reach of 
dark matter detectors, can be obtained in the scenarios analysed in this Subsection. 
The NMSSM nature is evidenced in these examples by the compositions of the lightest 
neutralino (which is a singlino-Higgsino mixed state) and the scalar Higgs (which can 
be mostly singlet and as light as m^o > 20 GeV). 

4.2 /lA^ < and /lAx < > 0) 

This choice comprises the cases /x, —Ax, —A^ > and — ^a, —A^ < 0. 

We first address the possibihty — ^a, —^k > 0. When compared with the cases 
discussed in the previous Subsection, the occurrence of tachyons in the Higgs sector 
gives rise to stronger constraints in this case, both in the CP-even and CP-odd Higgses. 

For CP-even Higgses tachyons are now more likely to occur, due to the negative 
contributions in 33, induced by the terms proportional to /lAi^ and Ax/ /i- Similarly, 
Ml II receives a sizeable negative contribution from the term proportional to /xAa 
which is particularly dangerous for large values of tan/3. In the CP-odd sector, an 
analogous study of the mass matrix shows that tachyons are more restrictive for large 
values of A and small values of k. Actually, from the naive requirement Alp n > the 
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following constraint is obtained n > —XAx/fi. In fact, this ensures the positiveness of 
the denominator in condition (ii) derived from the minimization of the Higgs potential 
in Section 2.2, and gives a qualitative idea on the dependence of the tachyonic region 
on the parameters Ax and /x. When compared with the cases treated in the former 
Subsection, larger regions of the parameter space are now excluded. We found that 
tachyons in the CP-odd sector typically give rise to stronger constraints than those from 
CP-even sector, although the corresponding excluded regions practically coincide. 

The experimental constraints from the neutralino sector are not very stringent in 
these examples, owing to the fact that the regions where the neutralino would have a 
small mass are typically excluded by the occurrence of tachyons. 

As an example, the (A, k) plane is represented in Fig. 11 for tan/3 = 3, 
Ax = —200 GeV, A^ = —50 GeV and /i = 110 GeV. In this case, and in contrast 
with what was displayed in Fig. 2, there exists a very large region where one cannot 
find minima 

of Wtfai- In particular, A > 0.25 is now excluded for this reason. In the 
rest of the parameter space, experimental constraints become very important in those 
regions with small values of the CP-even and CP-odd masses. Although the most im- 
portant exclusion is due to DHDM constraints {h^ hh and — > r+r^), some regions 
not fulfilling the bounds on APM {h^a^ — > 46's) also appear. It is worth emphasizing 
that in the remaining allowed regions the lightest Higgs is doublet-like (5^3 < 0.003) 
and its mass is never too small, m^jO > 85 CeV. 

Regarding the composition of the lightest neutralino, it turns out to be Higgsino-like 
in all the allowed parameter space {Nl^ + Nf^ > 0.98). As we already mentioned, those 
regions with small A and k that would lead to a singlino-like neutralino are excluded 
by the absence of physical minima in the potential. For this reason the mass of the 
neutralino is dictated by the value of the jj, term and we found m^o pa jj, throughout 
the allowed parameter space. 

This is shown in Fig. 11, which represents the corresponding values of the neutralino- 
nucleon cross section as a function of the neutralino mass. The cross section ranges 
from 6 X 10^^ pb < o'^o_p < 7 x 10^^ pb in this case. Once more, although higher 
values could be obtained, these are typically excluded due to the constraints on the 
CP-even Higgses. 

Let us now address the relevance of variations in A,^ in the allowed regions of the 
parameter space and thus on the predictions for (J^o^p. It can be seen that the increase 
in \Ak\ (i.e., making it more negative) translates into an almost negligible enlargement 
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Figure 11: The same as in Fig. 4 but for tan/3 = 3, = -200 GeV, = -50 GeV 
and fjL — 110 GeV. Only the lines with m^o = 114, 75 GeV are represented, and none of 
the lines showing the lightest scalar Higgs and neutralino composition is depicted, since 
S^s < 0.1 and Nf^ < 0.1 in all the plane. 

in the allowed area, while the experimental constraints on CP-even Higgses become 
more restrictive. On the other hand, a decrease in the value of l^^l leads to a lighter 
CP-odd Higgs and the tachyonic region increases, as can be easily understood from the 
mass matrix (2.10). For instance, in the particular case of — 0, and unless \Ax\ is 
also very small, the entire parameter space can be excluded. 

Regarding changes in Ax and /x, these clearly affect the regions excluded by tachyons. 
Large values of /x and small \Ax\ allow an increase in the accepted regions, in agreement 
with the condition on k derived above, k, > —XAx/fi. Also, note that, since the masses 
of the Higgses increase, the associated experimental constraints become less restrictive 
and the allowed area is larger. Nevertheless, the region where the neutralino would 
have an important singlino composition is still typically excluded, and therefore in the 
allowed region Xi is still Higgsino-like, with m^o ^ ^. Despite the increase of m^o, the 
predictions for the cross section are essentially unaltered. An example with tan/9 = 3, 
Ax = -200 GeV, A^ = -50 GeV and fi = 200 GeV is represented in Fig. 12, display- 
ing both the (A, k,) plane and the neutralino-nucleon cross section versus the neutralino 
mass. We find o-^o_p < 5 x 10~^ pb, similar to what was found in Fig. 11, but now 
with m^o pa 190 GeV. The singlino component of Xi is negligible (iVfg < 0.006) and 
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Figure 12: The same as in Fig. 4 but for tan/3 = 3, = -200 GeV, = -50 GeV 
and fj, — 200 GeV. Only the line with m^o = 114 GeV is represented, and none of the lines 
showing the lightest scalar Higgs and neutralino composition is depicted, since S^^ < 0.1 
and Nf^ < 0.1 in all the plane. 

the scalar Higgs is doublet-like (S'^g < 0.001). 

Finally, regarding variations in the value of tan /3, these have httle effect on the 
shape of the tachyonic region, whereas experimental constraints are more sensitive to 
them. As in the former scenario, for low values of tan P light scalar Higgses are obtained. 
Since these are predominantly doublet-like, experimental constraints (especially those 
associated with DHDM, namely, hP — > bb and hP — > r+r") become very important and 
forbid, for instance, the whole parameter space in the case tan/3 = 2. On the other 
hand, larger values of tan (3 are welcome in order to increase the value of m^^o , obtaining 
also a moderate enhancement of the cross section. In order to illustrate this discussion, 
we represent in Fig. 13 two cases with tan (5 — 2,5, for Ax — —200 GeV, = —50 GeV 
and yU = 110 GeV. We find that the cross section can reach cr^o_p 10^^ pb in the case 
where tan (3 = 5. None of these examples displays any qualitative change regarding the 
neutralino and Higgs compositions. 

To complete the analysis of the cases with ^Ax < we still have to consider the 
case fj,, —Ax, —A^ < 0. As we explained in the former Subsection, the analysis of the 
Higgs sector will be analogous to that of the case we have just studied. Despite the 
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Figure 13: The same as in Fig. 4 but for the cases Ax = —200 GeV, = —50 GeV, 
fj, = 110 GeV, and tan/? = 2, 5, from top to bottom. In the left frames, only the line with 
m/jO = 114 GeV is represented, and none of the lines showing the lightest scalar Higgs and 
neutralino composition is depicted, since S^^ < 0.1 and A^fg < 0.1 in all the plane. 
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differences in the neutralino sector, its mass, composition and detection cross section 
will also be qualitatively equal to those previously discussed. 

Summarizing, all the cases we have analysed in this Subsection present as common 
features the appearance of Higgsino-like ncutralinos with a detection cross section which 
can be as large as cr^o_p < 10~^ pb, and doublet-like Higgses. 

4.3 iiA^ > and /lAx > (k > 0) 

We consider now those cases where fiA^ > 0, and fiAx > 0, conditions which are 
fulfilled in the cases Ax, > and /x, ^4;^, ^4^ < 0. 

We will begin with all /i, Ax, and ^4^ positive. As in the previous cases, a simple 
analysis of the tree-level Higgs mass matrices gives a qualitative understanding on the 
nature and extension of the tachyonic regions in the parameter space. 

In this particular case tachyons in the CP-odd sector arise through the negative 
contribution — in Alp22- Since this is mainly compensated by the positive term 

\ 2 2 

^^Ax sin 2/3, the tachyonic region occurs for small values of A. The excluded region 
is obviously more important for small values of Ax and large /i, A^, and tan/3. The 
occurrence of tachyons in the CP-even Higgses is analogous to the case /jlAx > 0, 
IjlAi^ < discussed in Subsection 4.1, due to the increase of the off-diagonal terms in 
the mass matrix. As in that case, tachyons appear for large values of A and small k 
and become more stringent as tan /3 grows. 

Experimental constraints play also a very relevant role in this case. Close to the 
tachyonic regions the experimental bounds on the Higgs sector are very severe. In 
particular IHDM, DHDM {h^ — > bb, hP — > 2 jets) are responsible for the most impor- 
tant exclusions, although APM (mainly h^a^ — > 4 6's) may also be violated. Finally, 
excluding those regions where the direct neutralino production is in disagreement with 
the experimental bounds leads to important constraints in the region with light x?. 

An example with Ax = 200 GeV, = 110 GeV, A^ = 50 GeV and tan/3 = 3 is rep- 
resented in Fig. 14, depicting the constraints on the (A, k) plane and the corresponding 
predictions for (7j^o_j, versus the neutralino mass. In the small experimentally allowed 
region the lightest neutrahno is a mixed singlino-Higgsino state, with Nl^ ^0.4 and 
-^13 +-^14 ^ 0-6, and the lightest scalar Higgs can have an important singlet component 
{S^^ < 0.8). The experimental constraints impose m^^o > 70 GeV and m^jO > 85 GeV, 
which set a limit on the theoretical predictions for the neutralino-nucleon cross section 
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Figure 14: The same as in Fig. 4 but for Ax = 200 GeV, ^ = 110 GeV, = 50 GeV, and 
tan/3 = 3. Only the line with S^^ — 0.1 is represented, and none of the lines showing the 
neutralino composition is depicted since 0.5 > N^^ > 0.1 in all the plane. 

at cr^o-p < 4 X 10~^ pb. 

Variations of A,^ have an important impact on the allowed parameter space. As 
already commented, the region excluded due to tachyons in the CP-odd Higgs sector 
increases for larger values of For instance, in the example with Ax — 200 GeV, fi — 
110 GeV, and tan/3 = 3, the allowed region completely disappears for A^ > 110 GeV. 
On the other hand, decreasing the value of A^ the parameter space is enlarged. Recall 
that the minimal value ^4^ = has already been analysed in Section 4.1 in the context 
of a scenario with /i, Ax, —A^ > 0. 

Decreasing the value of Ax also leads to an increase of regions with a tachyonic 
pseudoscalar. li n — 110 GeV, A^ — 50 and tan/3 = 3, the whole parameter space 
is excluded for Ax < 50 GeV. On the other hand, a moderate increase of Ax helps 
avoiding tachyons, especially in the CP-even sector. An example with Ax — 300 GeV 
can be found in Fig. 15, where the (A, k) plane and the theoretical predictions for cr^o_p 
are represented. Since in this case the experimental constraints from Higgs decays are 
less severe, we find that the regions with very light Higgs and Xi a-re now experimentally 
viable. In particular, neutralinos with an important singlino composition, N^^ < 0.45, 
can be obtained with m^o > 45 GeV, whereas the hghtest Higgses (m^o fa 65 — 90 GeV) 
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Figure 15: The same as in Fig. 4 but for Ax = 300 GeV, ^ = 110 GeV, = 50 GeV, 
and tan/3 = 3. Only the lines with m^o = 114, 75 GeV are represented. Regarding the 
neutralino composition, only the line with Nf^ — 0.1 is shown in the upper corner of the 
allowed region, since in the rest of the parameter space 0.5 > N^^ > 0.1. 

are all singlet-like. This in turn favours larger values of the cross section {o'^o_p < 2 x 
10~^ pb), and compatibility with present experiments is almost obtained. Should we 
further increase the value of ^4^, the experimental constraints associated with the scalar 
Higgs would become again more important. 

In order to prevent the occurrence of tachyons in the CP-odd Higgses, the value of fi 
has to be small. For instance, taking = 200 GeV in the example with Ax — 200 GeV, 
A„ = 50 GeV and tan P — 3, all the (A, k) plane would be excluded. 

Regarding the value of tan (3, as already mentioned, the larger it is, the more ex- 
tensive the regions excluded by m^o < become. In Fig. 16 we represent two examples 
with Ax = 200 GeV, A^ = 50 GeV, = 110 GeV and tan/3 = 2, 5. We find that 
low values of tan /3 still allow physical minima of the potential. Moreover, both the 
singlino component of xl and the singlet component of the scalar Higgs can be slightly 
enhanced. For example, for tan P — 2 hght neutralinos (x? > 60 GeV) can be obtained 
with A^fg < 0.55, while Higgses in the mass range m/jO Ki 75 — 100 GeV are singlet-like. 
However, the predictions for the detection cross section suffer a moderate decrease and 
only cr^o_p < 5 x 10~^ pb is obtained. On the other hand, for tan /? > 5 the whole 
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Figure 16: The same as in Fig. 4 but for the cases A\ = 200 GeV, /i = 110 GeV, = 
50 GeV, and tan (3 = 2, 5, from top to bottom. In the case with tan (3 — 2 only the 
lines with m^o = 114, 75 GeV, and those with S^^ = 0.1 and Nf^ = 0.1 are represented. 
Similarly, in the case with tan/? = 5 only the line with nif^o = 25 GeV is drawn, and none of 
the lines showing the compositions of the lightest scalar Higgs and neutralino is depicted, 
since 0.7 > S^^ > 0.1 and 0.4 > N^^ > 0.1 in all the plane. 

parameter space is in general excluded. 

We should now address the complementary choice of the sign of the parameters, 
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namely /i, Ax, < 0, for which we aheady know that the analysis of the Higgs sector 
still holds. Once more, differences arise in the theoretical predictions for cr^o_p- 

To sum up, the choices for the signs of the parameters which have been considered 
in this section also permit obtaining large values for the theoretical prediction of the 
neutralino-nuclcon cross section, despite the fact that the parameter space is very 
constrained both experimentally and by the occurrence of tachyons. In particular, 
values of cF^o-p close to the sensitivities of the present detectors can be found in some 
regions of the parameter space. The lightest ncutralino displays a mixed singlino- 
Higgsino character, and the scalar Higgs is singlet-like and light in those regions with 
larger cr^o-p- 

4.4 iiA^ > and fiAx > {k < 0) 

We will now focus our attention on those cases with a negative value for k, namely 
/X, ^A, > and /X, ^A, < 0. 

Let us therefore concentrate on the first of the two possibilities, fi, Ax, A^ > 0. 
The parameter space is in this case plagued with tachyons in both the CP-even and 
CP-odd Higgs sectors. On the one hand, regarding the CP-odd Higgses, large values 
of \k,\ and fi and small values of A and Ax may lead to negative values in the diagonal 
terms of the mass matrix, especially in Ai^p^. Moreover, large values of A can also 
induce very large off-diagonal terms A^pi2 if Ax is large and a negative eigenvalue can 
be obtained in that case. Similar arguments lead to analogous conclusions concerning 
tachyons in the scalar sector, being the region with small A the one facing the most 
severe restrictions. 

Note that experimental constraints will play a very important role in the vicinity 
of these regions. Although the largest exclusions typically arise from the bounds on 
IHDM and DHDM (mainly in — > bb, 2 jets), APM can also exclude some 

regions with a small m^o. All these become particularly restrictive in the low tan/3 
regime. In fact, in most of the cases with tan/9 < 3 all of the parameter space is 
excluded. Experimental constraints in the neutralino sector can also be very stringent, 
especially for small values of /i, where Xi is light and Higgsino-like. 

In the remaining allowed regions of the parameter space, the lightest CP-cvcn Higgs 
is mostly dominated by the doublet component. Concerning the lightest neutralino, it 
turns out to be Higgsino-like. Owing to this, the predictions for (Jj^o „ are very similar 
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to those obtained in Section 4.2. 

The above discussion can be illustrated with Fig. 17, where the (A, k,) plane and 
the predictions for o'^o_p are presented for a case with Ax = 450 GeV, n = 200 GeV, 
Ak = 50 GeV, and tan/9 = 5. Light singlino-like neutralinos can only be obtained 
in the very small area with \k\ < 0.06, where the scalar Higgs may be as light as 
rrif^o > 50 GeV with a large singlet component. In this particular region the predicted 
values for the detection cross section are not large, cr.^Q_p < 10^^ pb. In the rest of the 
parameter space Xi is Higgsino-like, which implies m^o ^ ^, and the lightest scalar 
Higgs is a doublet with m^^o > 112 GeV. Slightly higher values for the cross section are 
obtained, which are bounded by the experimental constraints on the scalar Higgs at 
o'xi-p ~ 6 X 10^® pb. This prediction can be slightly increased with larger values of 
tan /3. For instance, with tan /5 = 10 one finds o'^o_p < 2 x 10~^ pb. None of the above 
remarks concerning the masses and compositions of the lightest neutralino and scalar 
Higgs would change in this case. 

Variations in the rest of the parameters are very constrained due to the extensive 
tachyonic regions and the strong experimental bounds, especially those associated to 
the bounds on IHDM and DHDM. This is, for instance, what happens when the value 
of Ax decreases. The very narrow region for small \k,\ where singlino-like neutralinos 
can be obtained is usually ruled out and the only surviving areas are those featuring 
heavy Higgsino-like neutralinos and doublet-like scalar Higgses. For this reason, the 
predictions for the detection cross section are always similar to those presented in 
Fig. 17. As an example of some of the most favourable results for the cross section 
that can be obtained, we show in Fig. 18 a case with ^4^ = 200 GeV, jj, ~ 110 GeV, 

= 200 GeV, and tan /3 = 10, where the sensitivity of the CDMS Soudan experiment 
is almost reached. 

Finally, concerning the case fi, Ax, < 0, nothing changes in the analysis of 
the Higgs sector. Once again, the differences in the neutralino sector and the slight 
changes in the experimental constraints can induce variations in the predicted o'^o_p. 
Nevertheless, a similar global upper bound of cr^o-p < 10"'' pb is obtained. 

The examples analysed in this Subsection feature a lightest neutralino which is 
Higgsino-like in most of the parameter space, together with a doublet-like lightest scalar 
Higgs. The ncutralino-nuclcon cross section is bounded by experimental constraints 
on the Higgs sector at o'^o-p ^ 2 x 10^^ pb. Singlino-like neutralinos can only be 
obtained in extremely small regions of the parameter space and predict smaller cross 
sections, cr^o_„ < 2 x 10^^ pb. 
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Figure 17: The same as in Fig. 4 but for negative k, with ^4;^ = 450 GeV, ji — 200 GeV, 
= 50 GeV and tan /3 = 5. 




Figure 18: The same as in Fig. 17 but for Ax = 200 GeV, = 110 GeV, A^ = 200 GeV 

and tan/3 = 10. 
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4.5 Variations in the gaugino mass parameters 

In order to complete our analysis we must now address variations in the gaugino mass 
parameters. These clearly affect the neutrahno sector, altering both the mass and 
composition of the lightest neutralino. In the former analysis we always assumed the 
relation /i < Mi < M2, which lead to neutralinos with important Higgsino com- 
positions. We will now generalize our results for different hierarchies among these 
parameters. Namely, we will investigate the consequences of having Mi < /i, M2 or 

M2 < fjL, Ml. 

The Higgs sector is not so sensitive to variations in the gaugino masses, since 
they only enter through loop corrections. For this reason, all the analysis regarding 
tachyons remains qualitatively valid, and the experimental constraints associated to 
Higgses exclude similar areas in all these cases. However, it is important to note that 
increasing the gaugino masses, especially the gluino mass, generally implies also an 
increase in the mass of the lightest scalar Higgs and for this reason, a heavy gaugino 
spectrum would typically lead to low values of o'^o_p. 

Let us first vary the values of the gaugino masses preserving their GUT relation, 
but allowing also changes in the /i parameter so that the relation Mi < /i, M2 can 
be achieved, thus increasing the gaugino character of x^. For this purpose we choose 
the example with — 200 GeV, = —200 GeV, and tan/3 = 3 that was rep- 
resented in Fig. 4, where now three different values for the /i parameter are taken, 
fjL — 110, 200, 500 GeV. Regarding the gaugino masses, we consider variations in the 
Bino mass as 50 GeV < Mi < 500 GeV, and the GUT relation Mi = | M2 = ^ M3. 
The results are shown in Fig. 19, where only those points fulfilling all the constraints 
are represented. Since the experimental constraint on the chargino mass imposes a 
lower bound on M2, the value of the Bino mass is also constrained. For this reason, 
for low values of // the lightest neutrahno is still a singhno- Higgsino state {Nf^ < 0.3 
and Nf^ + Nf^ > 0.7). However, as the value of fj, increases, so does the gaugino com- 
position of Xi (note that both the Higgsino and singlino compositions decrease when 
increasing fi). For instance, in the example with // = 500 GeV neutralinos lighter than 
"^xi ~ '^'^^ GeV are Bino-like. 

The appearance of Bino-like neutralinos has as a consequence the enormous decrease 
in the neutralino-nucleon cross section. In particular, the Higgs mediated interaction is 
now negligible and detection would only take place through the squark mediated inter- 
action. In contrast to the huge predictions for cr<>o in the case of a singlino- Higgsino 
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Figure 19: Scatter plot of the scalar neutralino-nucleon cross section as a function of the 
neutralino mass for Ax = 200 GeV, = -200 GeV, tan/3 = 3, = 110, 200, 500 GeV 
from left to right and from top to bottom. The gaugino masses satisfy the GUT relation 
and the gaugino mass parameter is varied in the range 50 GeV < Mi < 500 GeV. Only 
those points fulfilling all the constraints are represented. 

neutralino (see, e.g., the plot with /i = 110 GeV in Fig. 19), Bino-like neutralinos would 
have o"^o_p < 10^^ pb, and thus would be beyond the sensitivities of even the largest 
projected dark matter detectors. 
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Figure 20: The same as Fig. 19 but for the relation Mi = 0.2 M2. 

Very light Bino-like neutralinos can also appear when the GUT relation is relaxed. 
This is due to the freedom to choose very small values of Mi, unrelated to the more con- 
strained (from the bound on the chargino mass) M2. Furthermore, Bino-like neutralinos 
are not subject to such strong bounds on direct neutralino production as Higgsino- 
like neutralinos. In order to exemplify this possibility, we have chosen the relation 
Ml = 0.2 M2 (the GUT relation among the Wino and gluino masses, M2 = fMs, is 
preserved) and applied it to the same examples as before. The results are represented 
in Fig. 20. Neutralinos as light as m^o > 25 GeV are attainable, although the pre- 
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Figure 21: The same as Fig. 19 but for the relation Mi = 2M2 = | M3. 

dieted values for tlie deteetion eross seetion lie beyond the sensitivities of all present 
and projected experiments. 

Finally, more general compositions for the neutralino can also be found. In particu- 
lar, Wino-hke neutrahnos are attainable if M2 < Mi, This is an interesting possibil- 
ity, since Wino-like neutralinos predict in general larger values for the cross section than 
Bino-like Xi- order to explore this possibility we have repeated the same scan in the 
parameter space as above, but now with the relation Mi = 2 M2 — 1/7 M3 for the gaug- 
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ino masses at the EW scale. The results are shown in Fig. 21 for jj, = 110, 200, 500 GeV. 
Obviously, the Wino component of Xi becomes more important for large values of 
H. Thus, for instance, whereas in the example with n = 200 GeV Xi is a mixed 
Higgsino-Wino state (with ^ 0.6), Wino-like neutralinos populate the example 
with /J, — 500 GeV. Because of the experimental constraint on the chargino mass, 
Wino-hke neutralinos cannot be obtained below m^o 100 GeV. Note that although 
the theoretical predictions for the detection cross section are larger than in the case of 
Bino-like neutralinos of Fig. 19, there is an important decrease with respect to those in 
the case of Higgsino-singlino neutrahnos. In particular, in this case o-^o_p < 2 x 10~^ pb. 

All the variations discussed in this Subsection imply a variation in the mass and 
composition of the lightest neutralino and scalar Higgs. We have seen how Bino- 
and Wino-like neutralinos are attainable by decreasing Mi and M2, respectively, and 
increasing jj,. The singhno-Higgsino character of the neutralino is therefore lost in 
these cases. Furthermore, the increase in // also leads to doublet-like heavier Higgses. 
For these reasons the Higgs mediated interaction is no longer effective and the theo- 
retical predictions for the neutralino-nucleon cross section can have a huge decrease. 
Therefore, the optimal situation is the one we have analysed in all the examples of 
Subsections. 4.1 to 4.4, where is close to its lower accepted value and smaller than 
the gaugino masses. 

4.6 Overview 

In the previous Subsections we have presented a separate analysis of the distinct regions 
of the Ax, Ak, and tan/3 parameter space. This kind of approach was useful in order 
to comprehend the implications of individual variations of these parameters. Still, 
in order to obtain a global view on the theoretical predictions for cr^o_p, and their 
compatibility with present and projected dark matter detectors, it is useful to conduct 
a more general survey of the parameter space. Such an analysis constitutes a good 
overview of the properties of the cases studied in the previous Subsections. 

Since we are interested in regions predicting large cr^o_p, according to the conclu- 
sions of Subsection 4.5 we will focus our attention on the case // = 110 GeV, with heavy 
gaugino masses. Mi — | M2 — 500 GeV. The rest of the input parameters are allowed 
to vary in the ranges, -600 GeV < < 600 GeV, -400 GeV < A^ < 400 GeV, 
and tan/3 = 2, 3, 4, 5, 10, with A, \k\ e [0.01,0.8]. 

The results of this scan are shown in Figs. 22 and 23, for positive and negative 
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Figure 22: The same as in Fig. 3 but for ji — 110 GeV and the rest of the parame- 
ters in the ranges -600 GeV < Ax < 600 GeV, -400 GeV < < 400 GeV, and 
tan/9 = 2, 3, 4, 5, 10, and with positive values of k. Only those points fulfilling all the 
constraints are represented. 

respectively, where the theoretical predictions for CTj^o-p are represented as a function 
of the lightest neutralino and scalar Higgs masses. Fig. 22 therefore summarizes the 
results of Sections 4.1, 4.2, and 4.3. Points with large predictions for cr^o_p are found. 
These correspond to very light singlet-like Higgses, with even m^jO > 15 GeV, which 
are more easily obtained for low values of tan/3 (tan/9 < 5). The points in this region 
correspond to those with ^A\ > 0, discussed in Sections 4.1 and 4.3, respectively. The 
lightest neutralino in those points has an important singlino composition, Nf^ < 0.6. 

Fig. 23 generalizes the analysis of Section 4.4. As we had already discussed there, 
the lightest scalar Higgs is heavier, m^o > 75 GeV, and although it can be singlet-like 
in some cases, the predictions for the neutralino-nucleon cross section are typically 
low. Regarding the neutralino, it is Higgsino-like in most of the parameter space and 
therefore m^o ?a /i. 
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Figure 23: The same as in Fig. 22 but for negative values of k. 

5 Conclusions 

We have performed a systematic analysis of the low-energy parameter space of the Next- 
to-Minimal Supersymmetric Standard Model (NMSSM), studying the implications for 
the direct detection of ncutralino dark matter. We have thus computed the theoretical 
predictions for the scalar neutralino-proton cross section, fT^o_p, and compared it with 
the sensitivities of present and projected dark matter experiments. In the computation 
we have taken into account all available experimental constraints from LEP on the 
parameter space. 

Wc have found that large values of cr^Q_p, even within the reach of present dark mat- 
ter detectors (see e.g. Fig. 22), can be obtained in regions of the parameter space. This 
is essentially due to the exchange of very light Higgses, nif^o < 70 GeV. The NMSSM na- 
ture is evidenced in this result, since such Higgses have a significant singlet composition, 
thus escaping detection and being in agreement with accelerator data. In fact, Higgses 
as light as 15 GeV can be obtained. The lightest neutralino in those regions exhibits a 
large singlino-Higgsino composition, and a mass in the range 50 < m^o < 100 GeV. 

Let us finally mention that our work can be considered as a first step towards 
a more complete analysis of the direct detection of ncutralino dark matter in the 
NMSSM. As discussed in Sections 2.5 and 4, other potentially important constraints 
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on the parameter space should be addressed in forthcoming pubhcations. This is the 
case e.g. of the rehc density, the 6 — > 57 branching ratio, and the current upper hmit 
on the decay Bg — > 11^ ii~ . 
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A Relevant NMSSM interaction vertices 

A.l Higgs- quark- quark Yukawa coupling 

Parameterising the interaction of neutral CP-even Higgs fields with quarks as 

C,,h = -Qi [C'ylPl + C'yrPr] Qi K , (A.26) 

where /i° denotes the physical (mass) Higgs eigenstates and i — 1,2 up- and down-type 
quarks, one has 

2Mw sm jj 2Mw cos p 

where we have omitted the quark generations and S is the unitary matrix that diago- 
nalises the scalar Higgs mass matrix, defined in Eq. (2.9). 
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A. 2 Neutralino-neutralino-Higgs interaction 

The interaction of scalar Higgs and neutralinos can be parametrised as 

where a = 1 — 3 refers to the Higgs mass eigenstate, a,/? = 1 — 5 denote the physical 
neutralino states, and the couplings are defined as 

= {-9 (K, - tan ewK,) {SaiN;, - Sa2N;,) + 

- 2V2KSa3K,N;, , (A.29) 

= (A.30) 

In the text, and since we have exclusively analysed interactions involving the lightest 
neutralino states (i.e. a — j3 = 1), we have simplified the above as C^l = C^j^ and 



(A.28) 



A. 3 Neutralino-squark-quark interaction 

In terms of the mass eigenstates, the Lagrangian reads 

= 9. [CI'^'Pl + C^/'Pr] Xlqf , (A.31) 

where i = 1, 2 denotes an up- or down-type quark and squark, X — 1,2 the squark mass 
eigenstates, and a — 1, . . . , 5 the neutralino states. Since we have neglected flavour 
violation in the squark sector, only LR mixing occurs, and squark physical and chiral 
eigenstates are related as 




q q 
V2I V22 




(A.32) 



One can also make the usual redeflnition rjli — 77I2 = cos^qf and 77^2 = ~V2i — sinOg 
Therefore, for the up sector, and again omitting quark and squark generation indices. 
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the coefficients C£^* are given by: 



^ tan 9wK^ sin 9^ + ^ A^a4 cos 



-V2g 



AT* T" 



Q 



tan OwK^ 



cos + 



2Mwsmp' 



Cf 1 = -V2g 



— tan 9wN^^ cos - — — — 
2 2Mvi^ sm p 



A?^^4 sin 6lii 



sin + 



2Mvi/sin/3" 



A^*4 cos 



(A.33) 
(A.34) 
(A.35) 
(A.36) 



In the above, Fq(u) denotes the hypercharge of the SU{2)l quark doublet (up-singlet) 
and the isospin of the ul field. The analogous for the down sector is trivially 
obtained by the appropriate replacements ( Yu Yd. T^, rriu ma, 6u ^j, 

sin f3 — > cos/9 and Na4 Nas). In this paper, and since only q — q — Xi interactions 
have been considered, we have always used C}^^^ = Cf ]^, i.e., setting o; = 1 in the 
above. 
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